Simple shear is not so simple 
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Abstract 

For homogeneous, isotropic, nonlinearly elastic materials, the form of the ho- 
mogeneous deformation consistent with the application of a Cauchy shear stress is 
derived here for both compressible and incompressible materials. It is shown that 
this deformation is not simple shear, in contrast to the situation in linear elastic- 
ity. Instead, it consists of a triaxial stretch superposed on a classical simple shear 
deformation, for which the amount of shear cannot be greater than 1. In other 
words, the faces of a cubic block cannot be slanted by an angle greater than 45° 
by the application of a pure shear stress alone. The results are illustrated for those 
materials for which the strain energy function does not depend on the principal 
second invariant of strain. For the case of a block deformed into a parallelepiped, 
the tractions on the inclined faces necessary to maintain the derived deformation 
are calculated. 
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1 Introduction 



Batra pQ showed that in the framework of isotropic nonlinear elasticity a simple tensile 
load produces a simple extension deformation if the empirical inequalities hold. Here we 
consider the related problem: if we apply a shear stress to a nonlinearly elastic block, 
what deformation is produced? To formulate this question mathematically, we start by 
letting x(X) denote the current position of a particle which was located at X in the refer- 
ence configuration. Denote the unit vectors associated with a fixed Cartesian coordinate 
system in the reference configuration by (Ei, E 2 , E3) and in the current configuration 
by (e 1; e 2 , e 3 ). Consider now a rectangular block of a non-linearly elastic material, with 
edges aligned with Ei, E 2 , and E 3 . Batra pQ was interested in obtaining deformations 
x(X) associated with a Cauchy stress tensor of the form 

a = T(e 1 ®e 1 ), (1.1) 

where T is a constant, whereas here the concern is determining the deformations associ- 
ated with a Cauchy stress of the form 

<r = S(e l ®e 2 + e 2 ®e 1 ), (1.2) 

where S is a constant (in order to ensure that the equations of equilibrium are satisfied). 



We call (1.2) a uniform shear stress [21 [3], although it is also denoted "pure shear" in the 
literature, e.g. El El IE] • 

The problem being considered is therefore the stress formulation of the problem of 
shear (the "most illuminating homogeneous static deformation" according to Truesdell 
[9]). The approach adopted here is the inverse of the classical strain formulation of 
the same problem, first proposed by Rivlin [10]. In that approach, the relative parallel 
motion of two opposite faces of a block is considered. Motivated by linear elasticity and 
his own physical insight, Rivlin proposed that the resulting deformation of the block can 
be described by 

x = X + KY, y = Y, z = Z, (1.3) 

where (X, Y, Z) and (x, y, z) are the Cartesian coordinates of a typical particle before 
and after deformation, respectively, and K is the amount of shear. Using his now clas- 
sical constitutive theory, Rivlin was able to calculate the stress distribution necessary 
to support this deformation. This strain formulation of the simple shear problem has 
since received much attention in the literature with excellent summaries to be found, for 
example, in Atkin and Fox [6J and Ogden |llj . 

Note that in the linear theory, these two formulations are equivalent since an in- 



finitesimal shear stress of the form (1.2) leads to an infinitesimal deformation of the form 
(1.3) and vice versa. As noted by Rivlin [ID] , this equivalence does not carry over to the 
framework of finite elasticity. We address this problem here. 

We recall that the dead-load traction boundary condition in finite elasticity is the 
prescription of the first Piola-Kirchhoff stress vector on the boundary of the undeformed 
body. Here we focus on the application of a Cauchy shear stress, because it gives a one-to- 
one mapping between a <deformation, traction> pair and a stress vector on a boundary. 
By contrast, the assignment of a Piola-Kirchhoff traction leads to a many-to-one mapping 
and has thus an ambiguous physical meaning [12] . 
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2 Invertibility of the stress-strain relation 

Let B denote the left Cauchy-Green strain tensor. A homogeneous, isotropic, hyperelastic 
material possesses a strain-energy density which may be written as W = W(Ii, h, h), 
where h,h, h are the three principal invariants of B, given by 

(2.1) 



(2.2) 



h = tr B, I 2 = \ [If - tr(B 2 )] , I 3 = det B. 
The general representation formula for the Cauchy stress tensor er is then 

<T = f3 I + (3 1 B + f3_ 1 B-\ 

where 
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It is assumed that the so-called empirical inequalities hold, 

A) < 0, Pi> 0, < 0. 

If W is independent of h, then fl-i = and it follows from (2.2) that 

Assume now that < 0. Then Johnson and Hoger [13] have shown that 
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These coefficients are all positive by virtue of the empirical inequalities (2.4), 

ip > 0, Vi > °> ^2 > 0. 

For incompressible materials, only isochoric deformations are admissible and ^3 
at all times. Then W = W(h, h) only, and the representation formula is given by 
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where p is the indeterminate Lagrange multiplier introduced by the constraint of incom- 
pressibility and 

dW dW 
A =2_, ^ = -2^. (2.10) 

The empirical inequalities in this case have the form 

Pi > 0, p-t < 0. (2.11) 



The inverse form of this relation is given by (2.5), (2.6) with 0q replaced by —p. 

This invertibility of the classical stress-strain relation for both compressible and in- 
compressible materials plays a central role in what follows. 



3 Determination of the strain 



The proof that a uni-axial tension load produces an equi-biaxial contraction in isotropic 
finite elasticity is due to Batra pLJ. An alternative proof is given here, based on the 
invertibility relations of the last section. First note that a stress distribution of the form 
Oii.e. ' 



TOO 





(3.1; 



leads to the equi-axial deformation 
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as expected, where, if < 0, 

B n = ^ + ^iT + V 2 T 2 , 

and, if = 0, 

B 11 = (T-p )/p 1 , B33 



B 



33 



(3-3) 



(3-4) 



In both cases, Bn — B33 is of the same sign as T, owing to (2.8) and (2.11), which 
shows that universally, a uni-axial tension (compression) leads to a lateral contraction 
(expansion) . 

Following the same line of reasoning, we see that a Cauchy shear stress distribution 
of the form (1.2), i.e. 

'0 S 0" 

S 


leads to a left Cauchy-Green strain tensor of the form 



(3.5) 
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where 

B n = ip + ip 2 S 2 , B 12 = i(> 1 S, B 33 = ip . (3.7) 
if /?_! < and, if p_ x = 0, 

B u = £? 33 = -Po/Pi, B 12 = S/fa. (3.8) 
This is not compatible with the left Cauchy-Green strain tensor associated with the 



simple shear deformation (1.3) which has the form 



B 



1 + K 2 K 6 
K 1 
1 



(3.9) 



The strain (|3.7|) combines triaxial stretch with simple shear, a type of deformation 

Wineman 



which has been investigated previously by Payne and Scott [14] . Varga 
and Gandhi [16], Rajagopal and Wineman jTTj , and Destrade and Ogden p2]. To make 
this connection clear, we first note that Bu > 0, B\ x — B\ 2 > 0, -B33 > 0, because the 



strain tensor is positive definite. It follows that (3.6) can be written as 

X\ A2 y/l - XjX 2 2 

AlVl-AfA^ 
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(3.10) 



where Ai = a/ (B 2 ± — Bf 2 )/B ni A2 = y/Bn, and A3 = y/B^. (Here, and henceforth, we 
assumed that S > 0. In the case S < 0, a minus sign must be inserted in front of the 



square root.) A simple check shows that B 
can be decomposed as 



FF where the deformation gradient F 
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(3.11; 



corresponding to a triaxial extension with principal stretch ratios Ai, A2 > Ai, A3, followed 
by a simple shear with amount of shear y/l — XfX^ 2 - Here we notice a surprising universal 
result: any isotropic, hyperelastic material cannot be shared by an amount of shear 
greater than 1, independently of the magnitude of the shear stress; this limit corresponds 
to a maximum shear angle of tan _1 (l) = 45°. 
The corresponding deformation is given by 



x = XxX + X^l- AfA^r, y = X 2 Y, z = X 3 Z, (3.12) 
for compressible materials. The corresponding incompressible form is given by 

z = A-l ^A2 ^Z. 



X 
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(3.13) 



The principal stretches of the deformation (3.12), fi±, /i 2 and ^3 (say) can be computed 
from the general expressions found in [TTl [TH] - Here they are connected to the A's through 
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For infinitesimal strains, //j = 1 + £j, say, and the amount of shear K = yL — X 2 X 2 , 
becomes if ~ (ei + £2) 1 ^ 2 - At the lowest order in if (order 1), we find that Ai ~ 1, 
A2 — 1, A2if — K , A3 ~ 1 and thus, that the deformations (3.12) and (3.13) coincide 



with the small strain simple shear deformation (1.3). 



Finally, we note that since we have an expression for the deformation gradient F in 



(|3.11|), it is a simple matter to compute the components of the first Piola-Kirchhoff stress 

T 



tensor P = (det F) T (F x ) corresponding to the shear stress (1.2). We find that 
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(3.15) 



4 Special materials 

A wide class of isotropic materials is modeled by a strain-energy density W which does 
not depend on I2, the second principal invariant of strain. This leads to = and the 



stress-strain inversion is given by (2.5). 

It follows that -B33 = Bu in (3.6) and, further, that A2 = A3 in (3.11). In other words, 
for materials such that dW/dI 2 = 0, a shear stress (1.2) produces a simple tension in 
the shear direction (with tensile extension Ai and lateral contraction A2), followed by a 
simple shear of amount yl — A^A^ 2 , yielding a deformation of the form 



x = X 1 X + X 2 \/l-X 2 1 X 2 2 Y, y = X 2 Y, z = X 2 Z. 



(4.1) 



This is a two-parameter deformation, entirely determined by Ai and A 2 . These quantities 
are related to the shear stress and the material response functions as follows, 
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(4.2) 



where A2 is real by the empirical inequalities. 

In the incompressible case (the so-called generalized neo-Hookean materials), W = 
W(Ii) only, where ii is the first principal invariant. Then we have a one-parameter 
deformation, 
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Here, p is the Lagrange multiplier introduced by the constraint of incompressibility. Note 
that the strain-stress relation is now 



B 



V 
2W 



2W> 



-cr. 



(4.6) 



This yields, in particular, that p = 2W'B^, which is positive by the empirical inequalities 
and the positive definitiveness of B. 

Notice also from (4.3) that, in order for the deformation to be well-defined, < Ai < 1. 
Therefore for all generalized neo-Hookean materials, a shear stress leads to a simple 
shear combined with a uni-axial compression in the direction of shear (and thus, a bi- 
axial expansion in the transverse plane). In other words, the gliding plane Y =const. is 
"lifted" by the application of the shear stress (1.2), as can be seen in Fig. 1 (and the 
plane of shear Z = const, is pushed "outwards". 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



Figure 1: Cross-section of a unit block made of any generalized neo-Hookean material, 
subjected to an increasing shear stress of the form S"( e i <g> e 2 + e 2 <S> ex), where the e^ are 
aligned with the coordinate axes. The resulting deformation is a simple shear combined 
with a uni-axial compression in the direction of shear (here, successively of 0%, 20%, and 
40% along ex). The face initially at X = cannot be inclined beyond the diagonal dotted 
line at 45°. 



5 Normal tractions 



Somewhat counterintuitively in view of the assumed stress distribution (1.2), if one de- 
forms a cuboid (with faces at X = ±A, Y = ±B, Z = ±C, say) into a parallelepiped, 
as is usual when visualising simple shear, normal as well as shear tractions have to be 



applied to the inclined faces in order to maintain the homogenous deformation (3.12). 
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It follows easily from (3.12) that the unit normal, n, in the current configuration to 
the face originally at X = A is given by 



n = ( ni , n 2 , 0) = ^ _ _ 2 (^1, - V 1 - A?A 2 - 2 , Oj . (5.1) 

Let s denote the tangent vector defined by 

a = (-n 2j 7i 1 ,0). (5.2) 
The normal and shear tractions on the inclined face, originally at X = A in the unde- 



formed configuration, that are therefore necessary to maintain the deformation (3.12) are 
given respectively by 

\/l - A 2 A~ 2 

n = n • an = 2Sn\n% = —2S 2 1 _ 2 , 

2 — A X A 2 

S{n\-nl)=S 1 ^ ? . (5.3) 



s • an 



It follows immediately from these that 

- 2S < n < 0, < s < S, (5.4) 



and so the necessary normal traction on the slanted face is always compressive. 

For generalised neo-] 
simplify to the following 



—1/2 

For generalised neo-Hookean materials, A2 = X± , and these tractions therefore 



Af ' S - S/b ~2-Xi 



h = n/S = -2 \ \ § = S / S = —L-. (5.5) 



These forms are independent of the specific form of the strain-energy function and are 
plotted in Figure 2. Note, in particular, the rapid decay of the normal traction once 
loading begins. This suggests that a normal compressive traction of the same order as 
the shear stress has to applied to maintain homogeneity of the deformation once shear 
loading begins. Application of the shear traction is not as important. 

There are two important geometrical idealisations where these tractions on the in- 
clined faces can be ignored. The first idealisation is that of a thin block, for which 
8 > A. Here applying a constant shear stress to the faces Y = ±B generates the 



deformation (3.12). The second important idealisation is that of a half-space, with the 
semi-infinite dimension in the Y direction. Application of a constant shear stress then 
at Y = (and the generation of an equilibriating shear stress at infinity) again yields 



(3.12) 



Finally we conclude by recalling that Lurie [1] shows that for the shear stress (1.2), 
the two in-plane principal stresses are aligned with the diagonals X ± Y = const, and 
have both magnitude S/2. 
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Figure 2: Plots of the normal and shear components of the traction on slanted face versus 
Ai. 

6 Concluding Remarks 

Although all textbooks on nonlinear elasticity explain that a simple shear cannot be 
produced by shear forces alone, very few allude to what deformation is produced by a 
pure shear stress. In fact, we found only three such mentions, two correct: "In an isotropic 
non-linear elastic medium the state of pure shear is not accompanied by pure shear strain" 
[I] and "Pure shear is a significantly more complicated deformation compared with simple 
shear" [5] and one incorrect: "A state of uniform shear stress leads to a uniform (or simple) 
shear deformation" [3J . In this note we derived the homogeneous deformation consistent 
with the application of a shear stress. 

We wish to emphasize that this is not a purely academic problem. Indeed, in a 
laboratory it is the applied force that is easily controlled, not the deformation. The 
displacement in simple shear may be controlled only in special devices and in special 
situations. As Beatty [IHj put it: "In practice, [..] it is not likely that a global simple 
shear deformation may be produced in any real material". Similarly, Brown explains in 
his handbook Physical Testing of Rubber [20] that "With single and double sandwich 
construction, there is a tendency for the supporting plates to move out of parallel under 
load" , a phenomenon which clearly corrupts the concept of simple shear. As a remedy, the 
British Standard ISO 1827:2007 [21] advocates the use of quadruple shear test devices. 
In those tests, however, the distance between the plates is allowed to vary freely, and 
again, the simple shear deformation is modified. 

Therefore, the determination of the deformation field corresponding to a simple shear- 
ing stress field, in the framework of isotropic nonlinear elasticity, should be useful in the 
designing experimental protocols for real materials. We have shown that a pure shearing 
stress field produces a deformation field which is not the classical simple shear defor- 
mation, but a simple shear deformation superimposed upon a triaxial stretch, and that 
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the corresponding angle of shear cannot be greater than 45°. For the special case of in- 
compressible generalized neo-Hookean materials, a one-parameter isochoric deformation 
is obtained. For general isotropic materials the corresponding deformation field is more 
complex and even more complex for anisotropic materials (see, for example, [22| 123]). 
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